Hybrid quantization of an inflationary universe 



(N 

o 

(N 

>>: 



> 
On 



in 
o 

(N 



X 



Mikel Fernandez-Mendez 1 !*^ Guillermo A. Mena Marugan 1 !] and Javier 01medo{[] 
1 Institute de Estructura de la Materia, IEM-CSIC, Serrano 121, 28006 Madrid, Spain 

We quantize to completion an inflationary universe with small inhomogeneities in the framework of 
loop quantum cosmology. The homogeneous setting consists of a massive scalar field propagating in a 
closed, homogeneous scenario. We provide a complete quantum description of the system employing 
loop quantization techniques. After introducing small inhomogeneities as scalar perturbations, we 
identify the true physical degrees of freedom by means of a partial gauge fixing, removing all the 
local degrees of freedom except the matter perturbations. We finally combine a Fock description 
for the inhomogeneities with the polymeric quantization of the homogeneous background, providing 
the quantum Hamiltonian constraint of the composed system. Its solutions are then completely 
characterized, owing to the suitable choice of quantum constraint, and the physical Hilbert space 
is constructed. Finally, we consider the analog description for an alternate gauge and, moreover, in 
terms of gauge-invariant quantities. In the deparametrized model, all these descriptions are unitarily 
equivalent at the quantum level. 

PACS numbers: 04.60.Pp, 04.60.Kz, 98.80.Qc 



I. INTRODUCTION 

Several approaches attempt to combine the current 
classical description of the gravitational interaction with 
the principles of the quantum theory. One of the 
most promising candidates is the so-called Loop Quan- 
tum Gravity (LQG) [Ij , a non-perturbative, background- 
independent, canonical quantization of General Relativ- 
ity. Its application to cosmological models is known 
as Loop Quantum Cosmology (LQC) [2]. All the sys- 
tems which have been quantized in this framework show 
the same remarkable property: the classical big-bang 
singularity is replaced by a quantum bounce. This 
is indeed the case of massless scalar fields propagat- 
ing in Fricdmann-Robcrtson- Walker (FRW) scenarios 
with fiat H, closed [3, or open topologies @, with 
cosmological constant jg], or even in spacetimes with 
anisotropies Q- Also inhomogeneous systems have been 
studied, both perturbatively, considering tensor Q , vec- 
tor [§], and scalar Qj| perturbations, and non- 
perturbatively, as in the case of Gowdy cosmologies in 
vacuo 112-14] and with matter content Consider- 
able progress has been reached in situations in which 
the inhomogeneities propagate as a field in an effective 
background, where quantum corrections in the geometry 
have been partly incorporated. The analytic and numer- 
ical analysis carried out in Ref. 114] (treating the field 
classically) showed that, even in the presence of non- 
perturbative inhomogeneities (and disregarding the pos- 
sibility of extreme fine-tuning conditions), the quantum 
bounce persists. However, a more fundamental descrip- 
tion can be achieved when the field and the background 
are both treated quantum mechanically, so that, in par- 
ticular, their interaction is not limited to conform to any 
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particular effective dynamics [Tg|. This issue has been 
already addressed in Refs. [T^T [HI HDL making use of 
the fact that Gowdy scenarios with T 3 spatial topol- 
ogy, after a suitable gauge fixing, can be interpreted 
as a scalar field propagating in a Bianchi I cosmology. 
A completely quantum description of those models was 
achieved through a hybrid quantization scheme, proposed 
initially in Ref. fl2|, which combines LQC quantization 
techniques for the homogeneous sector with a standard 
Fock representation for the inhomogeneous degrees of 
freedom. In this approach, after a partial gauge fixing, 
the spatial average of the Halmitonian constraint is im- 
posed quantum mechanically, and the physical Hilbert 
space is constructed out of its solutions. 

Despite the variety of cosmological settings described 
above, more attention needs to be drawn to inflationary 
scenarios, given their fundamental role in the physics of 
the Early Universe [l7[ • Among all these scenarios, the 
case of a massive scalar field propagating in an FRW 
spacetime is of special interest: it is one of the simplest 
inflationary models keeping most of the required aspects 
for the satisfactory understanding of our universe. Re- 
cently, it has been possible to prove -at the LQC effec- 
tive dynamics level- that almost all of its solutions pro- 
vide enough inflation (l8| . solving the fine-tuning prob- 
lem arising in General Relativity. However, such promis- 
ing results do not follow from any genuine quantum dy- 
namics, owing to the lack of a complete quantization of 
the system. On the other hand, inflationary universes 
provide a natural framework for the development of pri- 
mordial cosmological perturbations [19]. The origin of 
those inhomogeneities could be explained by the assump- 
tion that they stem from the early vacuum fluctuations 
of the inflaton field. Assuming they were initiall y sm all, 
one can invoke perturbation theory to treat them [20l[21j . 
As far as observations are involved, scalar perturbations 
are the most interesting ones, since they are those that 
leave stronger imprints, e.g., in the cosmic microwave 
background. Thus, they have been observed with high 
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precision by the Wilkinson Microwave Anisotropy Probe 
(WMAP) Though, even when this type of inhomo- 
geneities has been partially analyzed in LQC inflation- 
ary models 

mum. 

it is essential to confront the results 
achieved so far with those that could be obtained from a 
complete quantization of the perturbed system, without 
any restriction to a particular effective background. 

To approach this problem, in this article we will carry 
out a complete quantization of an inflationary cosmolog- 
ical model with small inhomogeneities. In a first step, we 
will consider a homogeneous massive scalar field propa- 
gating in a closed FRW scenario [23[ . A quantum Hamil- 
tonian constraint is provided and its solutions are char- 
acterized. Specifically, it is noted that they are com- 
pletely determined by their initial data on the minimum 
volume section allowed by the discretization of the ge- 
ometry. The physical Hilbert space can be constructed 
out of this space of initial data just by equipping it with 
a suitable inner product. Once the homogeneous setting 
is established, we will introduce small inhomogeneities 
around the homogeneous solutions by means of pertur- 
bation theory. We will carry out a gauge fixing (adopting 
the longitudinal gauge), and a canonical transformation 
which includes the scaling of the matter perturbation by 
the FRW scale factor. The aim of such a transformation 
is to achieve a formulation of the system in which the dy- 
namical behavior of the matter perturbation approaches, 
in the ultraviolet limit, the one of a scalar field propa- 
gating in a static spacetime with an effective mass owing 
to the interaction with the FRW background. Next, we 
will proceed to quantize the system employing a hybrid 
approach, i.e., we will apply a polymeric quantization 
to the homogeneous sector while adopting a Fock rep- 
resentation for the inhomogeneities. The problem that 
then arises concerns the selection of a particular Fock 
quantization, as we have at our disposal infinite inequiv- 
alent representations. Fortunately, the recent u niq ueness 
results for fields in non-stationary spacetimes [251 ] moti- 
vate us to pick up the family of Fock representations in 
which the vacuum state is invariant under the group of 
spatial isometries -namely, the SO(4) group- and where 
the classical dynamics of the field (obtained after de- 
parametrization) is implemented by a unitary quantum 
operator. It has been possible to prove [26| that this 
family is indeed a unique unitary equivalence class, as in 
the situations treated in Ref. [25|]. This result is closely 
tied to the previously mentioned canonical transforma- 
tion, since a different (non-trivial) scaling of the field 
prevents the unitary implementation of its corresponding 
dynamics. Finally, we will formally provide the quantum 
Hamiltonian constraint of the system and its solutions, 
and endow them with a Hilbert space structure, thus ob- 
taining the physical Hilbert space. 

In order to show the robustness of the treatment and 
avoid problems related with the choice of a specific gauge, 
we repeat the procedure for an alternate gauge, in which 
the spatial metric is homogeneous. Moreover, in both 
situations we relate the modes of the resulting canonical 



pair of fields with a set of Bardeen potentials [2lJ , which 
are gauge invariants of the model, and prove that the 
canonical transformations relating our fundamental fields 
with these potentials are unitarily implementable in the 
quantum theory. 

This article is organized as follows. In Sec. [II] we pro- 
vide a classical description and the complete quantiza- 
tion of the homogeneous system. The inhomogeneities 
are included at the classical level in Sec. Mil and a com- 
plete hybrid quantization of the system is carried out in 
Sec. IIVI for a particular gauge choice. In Sec. [V] we re- 
peat the same analysis for an alternate gauge fixing. A 
description in terms of Bardeen potentials can be found in 
Sec. lVIl A summary with the main conclusions is given in 
Sec. IVIII Finally, to make the paper self-contained, four 
appendixes have been included. 

II. THE HOMOGENEOUS AND ISOTROPIC 

SYSTEM 

In this section, we will deal with the homogeneous and 
isotropic setting, which consists of a massive scalar field 
propagating in a closed FRW spacetime. We will give 
a classical description of the corresponding constrained 
system and proceed to quantize the model to comple- 
tion, determining the physical Hilbert space. Concerning 
the technical aspects of the quantization, we will mainly 
follow Ref. 0. 

A. Classical system 

The FRW spacetime under consideration can be foli- 
ated in spatial sections isomorphic to three-spheres. On 
them, we choose as fiducial metric the standard round 
metric of unit radius (instead of radius do = 2 as in 
Ref. [J]), with fiducial volume 1$ = 2n 2 . 

The degrees of freedom of the matter sector -without 
taking into account constraints- are a massive scalar field 
cj) and its canonically conjugate momentum p^, while 
the geometry will be described in terms of an Ashtekar- 
Barbero connection A l a — c °w*//o and a densitized triad 
Ef = pVU °ef/l%, with i being the internal SU(2) index 
and a the spatial one, both running from I to 3. Here, 
f2 a f, = °oj l a °cu l b is the fiducial metric of the spatial sec- 
tions, and il, its determinant. The classical algebra is 
given by {4>,p,p} = 1 and {c,p} = 8irGj/3, where G and 
7 are the Newton constant and the Immirzi parameter, 
respectively. 

Owing to the homogeneity and isotropy, only the 
Hamiltonian constraint H = NqCq/(16ttG) remains, 
where A^ is the (homogeneous) lapse function, 

Co = -^P[(c - lof + l 2 ll\ + ^-(pl + m W) , 

7 (2-1) 

and V = \p\ 3/2 . 
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B. Quantum system 



In addition, 



In the quantization of the system, we choose a standard 
Schrodinger representation for the scalar field, whose 
kinematical Hilbert space is "HJJ^* = L 2 (K, dtj)) ■ For the 
geometry, we apply a polymeric quantization The 
basic variables in LQC are fluxes of densitized triads 
through square surfaces, enclosed by four geodesic edges, 
which are basically given by p, and holonomies of the con- 
nection along integral curves of the fiducial triads of fidu- 
cial length fj,lo. The gravitational part of the kinematical 
Hilbert space is = L 2 (R Bo hr, e^Bohr), where R Bo hr 

is the Bohr compactification of the real line, and d/iBohr 
is the natural Haar measure associated with it. Then, 
the kinematical Hilbert space of the whole system is the 
tensor product <8> U^ v . 

To promote the Hamiltonian constraint to a quantum 
operator, its gravitational part must be written in terms 
of holonomies of the connection A l a . We will essentially 
follow the procedure of Ref . [j] , where the curvature ten- 
sor constructed out of su(2)-connections is regularized. 
As a first step, one constructs holonomies along the edges 
of a closed square, properly selected by considering four 
integral curves along alternating left and right invariant 
vector fields (left and right invariant vector fields com- 
mute), well adapted to the fiducial structures. Then, the 
curvature operator is basically replaced by a circuit of 
holonomies around a square enclosing a non-vanishing 
area A, determined by the infrared spectrum of the area 
operator defined in LQG (classically, the local curvature 
would be recovered in the limit A — > 0). As in Ref. 0], 
we adhere to the improved dynamics scheme, in which 
the fiducial length of the edges of the considered square is 
given by the function p, = y/ A/p, up to a factor Iq, In the 
triad representation, one can find a basis of normalizable 
states \v) in 'H5T V (with v £ K), on which the action of 
the matrix elements of the holonomies is Np\v) = \v+l), 
while p\v) — sgn(u)(27T7G'?iv / A|i'|) 2//3 |ti), with h being 
the Planck constant. 

Finally, we adopt the following operator as quantum 
Hamiltonian constraint: 



Co 



1/2 



8nG(p 2 + m 2 V 2 4> 2 ) - 



n 2 



^{(l+7 2 K^ 4/3 -^«in 2 (AW 



(2.2) 

1/2 



where V represents the volume V — |p| 3 ^ 2 , and we have 
introduced the inverse of this operator, defined as 



sga(v)V 

AirjGhvA 

1 3 

- NpV^ 3 N_p) 



(2.3) 




(2.4) 



ra(v) [e ^ N 2 pe » - e l » AT_ 2ri e' 



e 2 N 2 pe ^ 



e l i N_ 2 pe l ^ sgn(w) 



V 



1/2 



The constraint Co annihilates the state \v = 0) and 
leaves invariant its orthogonal complement. Therefore, 
we can decouple the state corresponding to the classical 
singularity, i.e., \v = 0), and restrict the study to its or- 
thogonal complement, that will be denoted from now on 
by V.^ v . Besides, with the usual definition of sin(/tZo) 
(see Ref. ||), fl 2 is the only operator in the constraint 
with a non-diagonal action on the i>-basis of "H|™ v . It 
only relates states with support on isolated points sep- 
arated by a constant step of 4 units in the label v and, 
moreover, different orientations of the triad are decou- 
pled. In conclusion, only states with support on semilat- 
tices £f — {v = ±(e + An); n £ N} are related by the 
action of ft 2 , where e £ (0,4] is a continuous parameter 
proportional to the minimum value of the physical vol- 
ume of the system in the sector under consideration. Let 
us emphasize that those sectors are preserved by all the 
operators of physical interest considered in this article. 
Consequently, they can be interpreted as superselection 
sectors, denoted from now on as 

In the following, and without loss of generality, we will 
restrict our study to % e + . Besides, we will apply a unitary 

transformation U = e doh ^ on where h(v) is defined 
in Appendix fAJ in complete parallelism with the analyses 
carried out in Ref. Q • This unitary transformation maps 
f2 2 into 



hi = un 2 u-\ 



(2.5) 



which is equivalent to the analog operator for a scenario 
with flat topology 



n -- 



i 



v 1 ' 2 



sgn(v)(N 2 p - N. 



4WA 

+ (N 2 p-N_ 2p )^(?)\v 1 / 2 , 



2/2 ) 



(2.6) 



just like the operator suggested in Ref. [2711 . but without 
inverse volume corrections (see Ref. [28] f° r additional 
information) . 

Hence, the gravitational part of the quantum Hamil- 
tonian constraint is finally encoded by 



1/2 





' l " 




' 1 " 


Cgrav 


v_ 


c 


V 



-1/2 



(2.7) 



with 



4- = -^|^-^sin 2 (^ ) + (l+7 2 K 2 ^ 4/3 }- (2- 
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As an aside, let us comment that, following the argu- 
ments of Ref. [J] (applied there to a related operator, 0), 
and noticing that the superselection sectors are now semi- 
lattices (instead of full lattices as in that work) , it seems 
reasonable to admit that the restriction of C gra v to H+ 
provides a self-adjoint operator with a non-degenerate 
discrete spectrum. 

The solutions (^\ annihilated by the constraint Cq can 
be defined in the dual of a dense set in <8> e.g., 
that of the product of the span of the w-basis and the 
functions of 4> with compact support. It is straightfor- 
ward to see that the coefficients &(v, </>) = (^\v, <p) satisfy 
a difference equation that resembles an evolution equa- 
tion with v playing the role of an internal time. Besides, 
one can check that the initial data ^(e,<fi) determines 
the whole solution at any volume v = e + An, Vn € N, 
in the considered semilattice. Therefore, we can iden- 
tify these initial data with the solutions, and complete 
the space of such data with an adequate inner prod- 
uct in order to reach the physical Hilbert space (29j . 
In this way, the physical Hilbert space can be taken as 
U phys = L 2 (R,d(f>). 



III. INHOMOGENEITIES: CLASSICAL 
DESCRIPTION 



In this section we will introduce small inhomogeneities 
in the model, treating them as perturbations. Specifi- 
cally, we will expand the Hamiltonian constraint up to 
second order in those perturbations. We will follow the 
same classical approach as in Ref. [3l[. Besides, we will 
carry out a gauge fixing in order to remove the unphysical 
degrees of freedom. Finally, we will introduce a suitable 
canonical transformation to prepare the system for its 
hybrid quantization, performed in Sec. IIVI 



Perturbations around classical homogeneous 
solutions. 



We now briefly revisit the inclusion of inhomogeneities 
in our cosmological model as perturbations of the homo- 
geneous system. For convenience, we introduce a per- 
turbative parameter e, making proportional to it each 
of the inhomogeneous corrections to the geometry and 
the matter field. Using this parameter, we carry out a 
perturbative expansion in the action, expressed in Hamil- 
tonian form, and truncate the expansion at second order. 
Furthermore, for simplicity, we will only consider scalar 
perturbations, something which is consistent at the stud- 
ied perturbative order because they are decoupled from 
genuine vector and tensor perturbations of the system in 
this approximation |3l| . Our focus on scalar perturba- 
tions is mainly motivated by their relevant role in obser- 



vational cosmology. In addition, we adopt a natural ex- 
pansion of the inhomogeneities in terms of modes of the 
Laplace-Beltrami operator in the three-sphere, namely, 
the (hyper-) spherical harmonics in S 3 : Q n , P™, and P™ b 
(see Appendix |B1). where we recall that Latin indices from 
the beginning of the alphabet denote spatial indices, and 
n € N + is here [HJ a positive integer that labels the eigen- 
values of the Laplace-Beltrami operator [see Eq. (|F3 ![) ] . 
Let us then write 



TV = 



aN (l + V2TT6 J^gnQ"), 



N a = cr 2 e a V27TeY,knPa, 



i 2 2 

h ab = a e 



n ab + 2 3 / 2 neJ2(a n Q n ^ab 



3« 



4> 



^fnQ n 



(3.1) 

(3.2) 

(3.3) 
(3.4) 



In these formulas, a 2 — 4itG/(31q), and the unperturbed 
scalar field has mass m = rh/a. In order to facilitate 
comparison with the anal ysis of Ref. [3l|, as well as to 
apply the results of Ref. [26| , we employ here the vari- 
ables a, if, and their corresponding momenta ir a and 
ir v , to describe the homogeneous sector. In the absence 
of inhomogeneities (i.e., in the unperturbed system), the 
relation with the homogeneous variables used in Sec. Ill Al 
is given by 



bl 



l 2 rr 2 p 2a 

<P 



p(c - l ) = -jl (T IT, 



,3/2 ' 



,3/2 
= t Q (JTV V . 



(3.5) 
(3.6) 



Note that the correspondence between a£l and the flux 
variable p is one-to-one, e.g., in the union of all the su- 
perselection sectors H+, to which we are particularizing 
our discussion. 

For the inhomogeneities, on the other hand, the canon- 
ically conjugate momenta of the coefficients a n , b n , and 
f n will be denoted by 7r are , ir bn , and irf n , respectively, 
while g n and k n play the role of Lagrange multipliers 
associated with two linear constraints. 

Truncating the action at quadratic order in e, we ob- 
tain the total Hamiltonian: 



H = iV 



H +e 2 Y,{H2+g n H$) + e 2 ^fe n if« I (3.7) 



where Hq corresponds to the Hamiltonian constraint of 
the homogeneous sector, 



-3a 



H 



I _2 | 2 i 6a ~ 2 2 4a\ f 

{-n a +TT +e m <p -e ), (3. 



and for each mode n we have 
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e -3a 
2 

(2a„7r a 



1 o 7l2 ~4 2 \ 2 / 15 2 n2 ~4, 2 \ 2 2 — 1 2 , 

-a n + 10-2— y6„ I ^ + I Y a « + 6 ^2T_~7 5 « J ^ - ""an + 4^- + ^fn ~ 6a n^f n ^ 

1 1 n 2 — 4 2 

-(n 2 - |) a 2 + - (n 2 - l)— -b 2 n + - (n 2 - 4)a„6„ - (n 2 - l)f, 



+ e bQ m 2 - 6^—-^ + ft + 6^a„/ f 



- 1 



iff 



iT\ = 



T? — 4 \ 

7T a — 7T a „ + 7T6 n + 3/„7T v 



K 7r <p 7r fn ~ 7r a 7r aJ _ (^q + 37r 2 )a„ - ^e 4a {(n 2 + l) " + _ 4 ) 6„ | + e 6a fh 2 (p(2f n + 3<pa n ) 



(3.9) 
,(3.10) 
(3.11) 



Here, £f 2 an d i^™ are the scalar constraints quadratic 
and linear in the perturbations, respectively, and iFj is 
the diffeomorphism constraint, also linear in the inho- 
mogeneities. Finally, note that, substituting N = <jN , 
the homogeneous constraint (|2.ip reproduces WttGHq as 
given in Eq. (|3.8[) . under the changes (|3.5I) and (|3.6I) . 



B. Gauge fixing 

For each mode n — 3, 4, 5, there are three dynami- 
cal degrees of freedom in the inhomogeneous sector: a n , 
b n , and /„, and two linear constraints: if ft and H^. 
The modes n = 1 and n = 2, on the other hand, will be 
treated independently, since they are in fact pure gauge. 
Hence, for n > 3, two non-physical degrees of freedom 
can be removed by means of a gauge fixing. Let us con- 
sider the longitudinal gauge, which is fixed by imposing 
the conditions 



o, n„ 



3tt v /„ = 0. (3.12) 



These conditions provide an admissible gauge since they 
are second class with the linear momentum and Hamil- 
tonian constraints. Besides, we can find values of k n and 
g n for which these gauge-fixing conditions are stable in 
the evolution. A simple computation yields that, on the 
gauge-fixed section, 



{b rll M} = 0^k n 
{n„,H} = 0^.g„ 



1 



n 2 - 4 



N e- 2a ir bn , (3.13) 
(3.14) 



Finally, if the conditions HJ\ = H\ = are satisfied, 
we can eliminate 7T6 ra , a n , and 7r a?i in favor of /„ and 
7Tf n . At this point we want to emphasize that, unlike the 
vanishing modes (b n ,TTi, n ), neither a n nor 7r 0re is equal 
to zero. Therefore, the reduced symplectic structure in 
the gauge-fixed system has a canonical form only in a 
suitable new set of variables on the reduced phase space. 



r 



Let us introduce the following set: 



a = a + e 



9 • 



ip = ip + e 2 ^3a„/„, 



fn = fn: 7f/„ = n fn ~ ^a n TT v , 

where, after reduction, it is understood that 



a„ = 3- 



e 45 (n 2 - 4) 



(3.15) 

(3.16) 

(3.17) 
(3.18) 

(3.19) 
(3.20) 



at the considered perturbative order. For simplicity in 
the notation, we have used a n in most of the expressions 
above, but one must keep in mind that its value is given 
by Eq. 

The reduced Hamiltonian constraint in this new set of 
variables, truncated to the correct perturbative order, is 



= A H + e- 



(3.21) 



where Hq has formally the same expression as Ho, but 
now in terms of the new variables a, (p, n a , and tt v [see 
Eq. (JHU)], and 



-3a 



EZL = 1 



nlE^ + f^E^ + f 2 ^ 



ff h 



3tt 



^2 



e 45 (n 2 -4)' 
2>TT a Ti v — e 6a m 2 (p 



(3.22) 



E 



ff 



f e 45 (n 2 -4) 
e 45 (n 2 -l) + e 6 "™ 2 -9^ 

3 



^2 



,(e 65 m 2 ^- 3n a n v ) 2 



e 4& (n 2 - 4) 
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Finally, we deal with the modes n = 1 and 2. For 
each of these cases, we only have two configuration vari- 
ables, a n and /„, and their corresponding momenta. 
Owing to the presence of the diffeomorphism constraint 
and the linear scalar constraint, these modes are com- 
pletely constrained. A convenient gauge fixing for them 
is a n = f n = 0. Then, the constraints HJ\ = Hj = 
imply that their corresponding momenta vanish. These 
conditions are stable under the dynamics if g n = k n = 0. 



C. Canonical transformation 

Now that the system has been reduced, we will change 
variables on the canonical phase space to adapt it to 
the requirements of the uniqueness results provided in 
Refs. [H, [2(|, regarding the quantization of fields in 
non-stationary scenarios after deparametrization. Specif- 
ically, any S'0(4)-invariant Fock representation that im- 
plements the dynamics unitarily is only compatible with 
a particular choice of variables on the phase space of the 
field, as any genuinely time-dependent linear canonical 
transformation of those variables prevents the simultane- 
ous fulfillment of both properties in another representa- 
tion [33[ . We can introduce that preferred set of variables 
by means of a canonical transformation on the reduced 
phase space. Let us start with: i) a scaling of the field 
of modes /„ by the FRW scale factor, and ii) the inverse 
scaling of its momentum, also allowing a suitable mo- 
mentum shift proportional to the configuration variable. 
These modifications can be extended straightforwardly to 
a canonical transformation which involves both the ho- 
mogeneous and the inhomogeneous sectors. In this way, 
we arrive at the following variables, which are canonical 
at the considered perturbative order: 



n 

fn = e & f n , w fn =e- 5 (7f /n -TrJ n ). (3.23) 

After the change, the homogeneous part of the Hamil- 
tonian constraint remains formally the same, though with 
the old homogeneous variables replaced with the new 
ones, while the perturbed Hamiltonian at order e 2 is 

m = ^ A + + /«%) > ( 3 - 24 ) 
EHL = 1 * 



e 4a (n 2 -4)' 



™ _ R - 2ir a ii lp - e 6a m 2 <p 
- bn v — 



e 6a (n 2 - 4) 



E? t = n 2 - 1 - 



vf 2 + 15tt 2 - e 4a + 3e 6a mV 



Jff " ' 2r'" 

2a - 2 o ( 2 ^ v - e m m 2 (p) 2 



This description of the reduced phase space, and of 
its Hamiltonian constraint, will be our starting point for 
the hybrid approach, except for a convenient transforma- 
tion in the homogeneous sector of the phase space, car- 
ried out to reintroduce a natural set of variables for its 
polymeric quantization. This transformation leads from 
a, it a , <f>, and n v to variables that are formally similar 
to p, c, 4>, and p^, namely the variables introduced in 
Sec. Ill Al though defined now in coexistence with the in- 
homogeneities. The change is given again by Eqs. (|3.5j) 
and p.6[) . replacing the old variables with their barred 
counterpart. Besides, we will call Co and , respec- 
tively, the constraints Ho and expressed in terms of 
these new variables for the homogeneous sector, and mul- 
tiplied by a factor of 16irG — 12ZqCT 2 in order to adopt 
the usual conventions employed in LQC. 



IV. HYBRID QUANTIZATION 

As we have anticipated, our main aim is to accomplish 
a complete quantization of this inflationary model con- 
taining inhomogeneities. As a preliminary step, we will 
introduce an auxiliary quantum framework -the so-called 
kinematical Hilbert space- to construct a quantum rep- 
resentation of the classical system, and in particular, of 
the Hamiltonian constraint. In order to determine the 
physical sector, we will look for the states annihilated by 
this quantum constraint, and build the physical Hilbert 
space out of them. 



A. Kinematical Hilbert space 



Let us consider again the kinematical Hilbert space 
H^n ® ^idrT f° r tne homogeneous sector, as explained 
in Sec. IIIBl We recall that, in %Jgf, the operator $ 
acts by multiplication and its conjugate variable p^ is a 
derivative operator. In the gravitational sector, the fun- 
damental variables are fluxes and holonomies of su(2)- 
conncctions, essentially represented, respectively, by the 
multiplicative operator p and the matrix elements of the 
holonomies, i.e., Np (in the improved dynamics scheme). 
Both operators have a well defined action on HST V . 

For the perturbations, in Sec. IIII CI we arrived at a 
privileged description, given by the variables f n and Wf n , 
defined in Eq. p.23[) . In this description, the massless 
representation permits a unitary quantum implementa- 
tion of the field dynamics (26|. To perform a standard 
Fock quantization, it is hence convenient to rewrite these 
modes in terms of the corresponding annihilation and 
creation-like variables, namely at n and their complex 
conjugate a*^ , with 



1 



-(oJnfn + WT/ n ), 



(4.1) 



e za m z - 3- 



e 8a( n 2 _ 4) 



(3.25) 



and w 2 = n 

operators &f n and Oj , such that 



1. We promote these variables to quantum 
fW rA. () at J = S nn >. Let 
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us now call 5 the vector space consisting of finite linear 
combinations of the A-particle states 



|A0 = \N 3 ,N A ,...), ^iV„<oo, 



(4.2) 



n>3 



with N n € N being the number of particles of the n-th 
mode. Then, we can construct the inhomogeneous sector 
of the kinematical Hilbert space, J 7 , as the completion 
of S with respect to the inner product (Af\Af') = #aW- 
Obviously, the A-particle states are then an orthonormal 
basis of the Fock space J- '. The total kinematical Hilbert 
space can thus be taken as 



u 



kin 



u 



kin 



n 



kin 



> T. 



(4.3) 



On this Hilbert space, it is understood that both df n and 
a , act as the identity on the homogeneous sector. 



B. Quantum constraint 

Let us now concentrate our discussion on the Hamil- 
tonian constraint Q3.24|) after rewriting it in terms of the 
variables p, c, <fi, and p^, as we have commented above. 
To introduce a(n at least symmetric) operator that rep- 
resents this constraint on H kln , we admit the natural as- 
sumption that the superselection sectors of the homo- 
geneous system (i.e., the sectors W+, if we are already 
restricting the analysis to positive p) are preserved by 
the inclusion of the inhomogeneities, since the latter are 
regarded as perturbations of the homogeneous setting. 
We then adopt the following quantization prescriptions: 

1. Contributions like <j)p<p will be promoted to the op- 
erator (fip^ +Pfj,4>)/2. 

2. In any product of a power of the volume V and 
a non-commuting expression, the former will be 
evenly distributed around the latter so as to ob- 
tain a symmetric combination. 

3. Any even power of the form [(c— lo)p] 2k , with k £ Z, 



will be represented by O 



(fc) 



= [Q '?] , constructed 

using the positive operator fi 2 [see Eq. (|2.4j) ] and 
the spectral theorem to define its fc-th power [34], 

4. In the case of odd powers of the form [(c— lo)p] 2k+1 , 
the prescription will be 



l )p] 2k+1 -+ &° (k) = |fi| fe A|f^, (4.4) 



A = 



where |fi| = V Q 2 and 

V 1 ' 2 



8WA 



[sgn(u) { r 



r=+l,-l 

J2 (re~ lriil0 N 4rp ,e- ir ^ h ) Bgpjv)] V 1/2 . (4.5) 

r=+l,-l J 



One can straightforwardly see that A is a difference 
operator which only relates states with support in semi- 
lattices of step 4 for v. Therefore, any power of (c — la)p 
will be promoted indeed to an operator that preserves 
the superselection sectors in this homogeneous volume 
(see Sec. EH). 

Using these assignments, we arrive at a quantum con- 
straint of the form 



(4.6) 



where Cq is the homogeneous constraint defined in 
Eq. (|23), and 



c 2 n = eiy 



1 

2uj„ 



V 6 r- 



" 1 " 












y _ 










3iw, 






f — 



1 



-X1G 



(4.7) 



1/6 




At this stage, it is worth noticing that any possible factor 
ordering ambiguity affecting C^, resulting from prescrip- 
tions other than ours for its operator representation, pro- 
duces only subleading quantum geometry corrections to 
the (already) perturbative terms in the total constraint. 

As in the absence of inhomogeneities, the constraint 
C annihilates the state \v = 0) (times any state in 
^mat ® -T 7 )' an d leaves invariant its orthogonal com- 
plement, 'H kln , complement to which we can hence re- 
strict all considerations in the following, removing the 
state analog of the cosmological (homogeneous) singu- 
larity. Besides, as anticipated, the operator C preserves 
the sectors of superselection in the homogeneous volume, 
namely the (positive) semilattices of step 4 in the label v, 
since all the basic operators from which C is constructed 
preserve those sectors in fact. In other words, as in the 
analysis of the homogeneous system, the action on v of 
all the operators considered in this article preserve the 
subspaces T-L e + of states with support in the semilattices 
Cf (limiting again the discussion to positive values of v 
for simplicity). 
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C. Physical Hilbert space 

In order to complete our quantization, we still have 
to characterize the space of solutions to the constraint 
C and provide it with the structure of a Hilbert space. 
In accordance with our perturbative approach, we will 
assume that the solutions can be expanded in a pertur- 
bative series and truncate them in the form 



m = w<i (0) 



(2) 



(4.11) 



In turn, each term can be expanded employing the basis 
that we have introduced for the total kinematical Hilbert 
space of the system: 

w<i (fc) = E E / ® (v\ ® (<^# (fc) (A^,</>). 

. (4-12) 

By consistency, the Hamiltonian constraint (ip\C must 
vanish order by order, up to the level of our approxima- 
tion in the perturbative expansion. Here, the dagger de- 
notes again the adjoint. The zcroth-order contribution 
to the constraint yields (V'l^C'd = 0. But this condi- 
tion is just the constraint studied in Sec. Ill Bl for the 
homogeneous sector. Recall that we have shown that the 
corresponding solutions are completely characterized by 
the initial data on the minimum volume section v = e. 
In addition, note that, at this zeroth order, the informa- 
tion that (ip]^ codifies about the inhomogeneities does 
not change in the evolution in v that the constraint Cq 
dictates. 

The next term in the perturbation expansion gives 



M( 2 )C t + (#o)(£^) + =0. 



(4.13) 



This relation tells us that ip^ must satisfy a difference 
equation similar to that for (ifjy-°\ but now with a source 
term which accounts for the interaction of the inhomo- 
geneities with the "background" state (ipl^ ■ Again, this 
can be interpreted as an evolution equation in the inter- 
nal time v, in which any solution emerges out of the ini- 
tial section without the need for any boundary condition 
around it. Consequently, if the initial data ip^\ v=E is 
provided, one can straightforwardly determine ip^ \v=e+4 
once (ip\( ' is supplied: one only needs to solve a (well- 
posed) linear difference equation in v. By the same ar- 
guments, following an iterative process we can find the 
value of the solution on any other section v — e + 4n, 
Vn G N. In conclusion, any solution to the constraint 
can be determined if the initial data on the minimum vol- 
ume section (v = e) are given. Although the construction 
of the solutions is formal, in general, we will see that this 
suffices to determine a physical Hilbert space which re- 
tains all the true degrees of freedom of the system. 

To do that, we will proceed in a similar way as we 
did for the homogeneous part of the model in Sec. Ill Bl 
Again, the important point is the fact that we can iden- 
tify solutions with their initial data on the minimum vol- 
ume section. Therefore, we just need to equip that space 



of initial data with a suitable inner product in order to 
construct the physical Hilbert space (see Ref. 13] for ad- 
ditional discussion). This inner product can be fixed, for 
instance, by demanding reality conditions on a complete 
set of observables [3(|. Implementing this approach, the 
physical Hilbert space that we obtain is then 



V. AN ALTERNATE GAUGE FIXING 

In addition to the previous analysis, we will consider 
now an alternate gauge fixing in which the spatial metric 
looks homogeneous (and hence also the spatial curva- 
ture). We will carry out the reduction and quantization 
of the system in a similar way to what we did in the 
gauge discussed so far. 

Let us start again with the constraint given in 
Eqs. (|3~TU|) . and (|3~TTj) . We will impose the condi- 

tions a n = and b n — 0, which are second class with the 
linear constraints. The stability of this gauge under the 
evolution can be ensured by fixing g n and k n as certain 
linear functions of the momenta n an and TTb n , whereas 
these momenta can be written as functions of the vari- 
ables /„ and 7r/ n -and the homogeneous variables- once 
the constraints f/]™ = and H\ = are satisfied exactly. 

Specifically, requiring that {b n , H} = 0, we arrive again 
at the very same condition given in Eq. (|3. 131) . while 
(together with this last equation) the consistency demand 
{a n ,H} = leads to 



1 

9n = 

7Try 



n z - 1 



-7Tfc„ - 7T Q „ 



(5.1) 



where we are using the following relations, obtained in 
the reduction process: 



tt v m 2 e 6a tp 

7Ta„ = 7T/„ H Jn 

7T„ 7T„ 



(5.2) 



7Tm 

n b n — T/„ 



m 2 e 6a ip 



After introducing these last expressions in Eq. p.9[) , we 
arrive at a Hamiltonian constraint of the same type as in 
Eq. (|3.2ip . but with the new coefficients 

3 7T 2 



7T„ W- 



TT a IT, 



iLfj = e [n — 1) + e m +9^ 

3 / e 6a mV x2 

37T W 



(5.3) 



It is evident from all these equations that the intro- 
duced gauge fixing is always well posed except on the 
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section of phase space corresponding to vanishing mo- 
mentum TT a . We will see later that even the potential 
problems posed on that section are eluded in our quanti- 
zation, owing to our regularization prescription and the 
fact that the kernel of the quantum counterpart of 7r a be- 
longs to the continuum spectrum, so that its correspond- 
ing operator can be inverted via the spectral theorem 

We will now introduce a scaling of the configuration 
variables f n by the FRW scale factor, extended to a com- 
plete canonical transformation, so that the new canonical 
pair of fields that describe the matter perturbation admit 
a Fock quantization with the good properties explained 
in Sec. IIII CI To this end, in addition to the scaling of 
the /n's, we perform the inverse scaling of the conjugate 
momenta, to which we also add a term linear in their cor- 
responding configuration field variables in order to ensure 
that the cross-term coefficients E^ have a subdominant 
contribution to the matter field dynamics in the large n 
limit. Of course, we must transform the homogeneous 
variables as well, so that the entire transformation on 
phase space is canonical at the considered perturbative 
order. Explicitly, the canonical change is given by 



<p = <p- 



fn — e a fn, 7T/„ 



1 cp " ip ) 

n 2 



3— + TT a ) fn 

7T Q 



(5.4) 



Under this transformation, the homogeneous part of 
the Hamiltonian constraint (|3.8[) is kept formally the 
same, though with the old homogeneous variables re- 
placed with their new counterparts, whereas, in terms 
of the new inhomogeneous variables, the contributions 
quadratic in the perturbations have again the general 
form given by Eq. (|3.24[) . but now with coefficients 



El 



1 



4tt 2 



E 



ft 



4) 



e Qa rh 2 (pir v 



E n f * = n z -I- — - 



+ e 2a m 2 



+ 



3e 2 



-rh 2 ip 



30tH 



, 3# 



- 1 



4) 



e ea rh 2 (p 



27f„ + 



(5.5) 



Finally, let us consider the variables c, p, <p, and p<p 
applied in the quantization of the homogeneous system, 
already described in Sec. [TTJ Again, they are given by 
Eqs. (|3.5p and (13.61) . replacing the old variables with 
their barred counterparts. Once more, we can combine 
a polymeric quantization of this homogeneous sector to- 
gether with a standard Fock quantization of the inho- 
mogeneities, providing a kinematical arena for the quan- 
tum treatment of the system. Following the quantization 
prescription of Sec. IIVB1 the reduced Hamiltonian con- 
straint, with coefficients given by Eq. (|5.5I) . can be pro- 
moted to an adequate operator, like the one introduced 
in Eq. (14.61) and with a similar contribution (|4.7|) of the 
inhomogeneities, but now with the operators 



Ft = - 



3a 2 l 



+ 



<J 2 k 



m 2 \ P \ 1/2 (-4 7 (# +W0)e?_ 1) - -7V/ 3 , 

2/3 / 



e 



Q„,2^.2;3 -3 Ae 
37 CT « O P0 B (-1) 



30p£ - 277^/^9^ 

2/3 

1 ? 

+ 2 + -l 



(1) 



2 llo 2 



7^6 







2/3 

3<J 2 l 




2/3 


) 


v _ 


+ 

n 2 - 4 


y 


( 



r 


■ i ■ 




V 



1 m 2 9 n 2 — 1 , 



'o 



\p\ i/2 + i%pl®U)\ 

2^ + 7 m 2 |p| 3 / 2 e^ 1) |p| 3 / 2 

(5.6) 



G = 2 7 Z 2 



1/3 



" 1 " 




y _ 





m+pMp\ 3/2 ®U)\p\ 3/2 + iffii-x) - ^%p%\p\ 1, ^U)\^' 2 



1/3 



This Hamiltonian constraint also decouples the zero vol- 
ume state (tensor product any state on the matter field, 
both for its homogeneous part and its inhomogeneities), 
and is a combination of operators that preserve the super- 
selection sectors of the homogeneous sector: semilattices 
of step 4 in the label v. On the other hand, notice the 



appearance of negative powers of Cl 2 through the terms 
©(_!)) ®(-2)' an d ®(-iV This inverse powers are well 
defined via the spectral decomposition of il 2 , since the 
kernel of this operator is in its continuum spectrum [35j . 

Physical states are annihilated by the operator C. We 
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will assume the same perturbative form for them as 
in Eq. (|4.11l) (at the order studied in our approxima- 
tion). The part of the solutions represents again 
a background-like state. Following the analysis carried 
out in Sec. [Ill anv such solution (at zeroth perturbative 
order) is characterized by its initial data on the min- 
imum volume section. Once that background solution 
is determined, one can evaluate (ip\^ 2 \ which takes into 
account the interplay between the inhomogeneities and 
the homogeneous background. In Sec. HVCl we saw that 
this part of the solution can be completely determined 
from its initial data at v = e once (i/j|(°) is known, at 
least formally. In total, we conclude that the space of 
solutions to the constraint C can be identified with the 
space of data on the minimum volume section. Finally, 
we can select a convenient inner product on this space, 
e.g., by choosing a complete set of observables and de- 
manding that the complex conjugation relations between 
their elements become adjointness relations between their 
corresponding operators. In this way, we reach the same 
conclusion as in Sec. lIVCl the physical Hilbert space can 
be taken unitarily equivalent to the Hilbert space (|4.14|) . 



VI. GAUGE-INVARIANT FORMULATION: A 
UNITARY MAP 

Gauge-invariant quantities provide a physically mean- 
ingful description of cosmological perturbations [2lj |. 
Such quantities are usually employed to describe the 
physics in a consistent manner, independent of the iden- 
tification of the spacetime and its matter content when 
transformations under diffeomorphisms are allowed, and 
insensitive to the specification of a particular gauge. In 
order to consolidate our proposal and show the robust- 
ness in this respect, in this section we will establish the 
correspondence between our fundamental variables /„ 
and itf n and gauge- invariant quantities, thus reformu- 
lating our quantum description in terms of the latter. 
Furthermore, we will see that both formulations (namely, 
the original one in terms of the matter field perturbations 
and the new one in terms of gauge invariants) are related 
by a canonical transformation which can be implemented 
quantum mechanically as a unitary transformation, at 
least, after deparametrization of the theory (e.g., in the 
regime of quantum field theory in a curved background 
spacetime). 

Let us consider a transformation of the general type 
e£ M , where x^ is a spacetime point, £ M is an 



On the other hand, instead of following the procedure 



arbitrary vector, and x' 1 * is the transformed point (again, 
e is the perturbative parameter introduced in Sec. MI Al . 
We express the covariant counterpart of £ M using (hy- 
per-) spherical harmonics, obtaining the mode expansion 



n 



(6.1) 
(6.2) 



presented in Ref. (11| to obtain gauge-invariant canon- 
ical pairs, here we will rather start from the standard 
formalism of Ref. (2l| , where several gauge- invariant per- 
turbations are defined. In particular, we will consider the 
gauge- invariant quantities $„, £™, and given by 
Eqs. (|C3 ]> -([C6 ]l of Appendix O The first couple of invari- 
ants is constructed purely out of geometrical quantities. 
The second one contains also matter perturbations, and 
we will refer to them as the energy density and the ve- 
locity gauge-invariant perturbations, respectively. Any 
other gauge-invariant variable can be expressed as a lin- 
ear combination of these four ones. This is the case of the 
Mukhanov-Sasaki variable (see Appendix ID]) . which is a 
gauge-invariant perturbation well adapted to flat scenar- 
ios, where it satisfies the equation of motion of a scalar 
field propagating in a static spacetime but with a time- 
dependent potential (see, e.g., Ref. [ll[ for a recent dis- 
cussion). However, for the matter content considered in 
this article, the gauge-invariant energy E™ (up to a suit- 
able scaling) is more convenient inasmuch as it satisfies 
the same type of equation but independently of the spa- 
tial topology, a fact that makes it privileged in compari- 
son with other invariants. 

In fact, one can check that the combination of 
Eqs. ([CTT]) and (jCT2l) with (|CT6| and (ICT71) yields a sys- 
tem of first order differential equations for the quantities 
£™ and v * which resembles the first order equations of 
a canonical pair of variables describing an oscillator with 
time-dependent frequency. Taking these considerations 
as our starting point, we define the following variables 



5a 



n*„ = - 





-4 


vVi 2 - 


- 4 TT 


\Jn 2 - 


- 1 e< 


e 3a 


1 




-4 



in s 

V 



e 6a rh 2 tp 



27T Q 



(6.3) 



Using Hamilton's equations for 7r v and a, one can ac- 
tually see that the two introduced variables are related 
dynamically by 



(6.4) 



where the dot stands for the derivative with respect to 
the conformal time (corresponding to No = e a ). With 
the above equation of motion and the corresponding one 
for , one obtains a second order differential equation 
for ^ n that corresponds exactly to that of the modes of 
a scalar field propagating in a static spacetime with a 
time-dependent quadratic potential. 

In the rest of this section, we will show that, for each of 
the two gauge fixations considered in our discussion, the 
pair of gauge invariants defined in Eq. (|6.3I) are related 
with the fundamental fields /„ and 7tt n by means of a 
canonical transformation. Moreover, we will also prove 
that such a canonical transformation is implementable 
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as a unitary transformation in the corresponding (de- 
parametrized) quantum theory. 



A. Longitudinal gauge b„ — H n = 

In Sec. IIIII we analyzed this gauge and introduced a 
canonical pair of variables /„ and Kf n , given in Eq. (|3.23|) . 
to describe the inhomogeneities after the corresponding 
reduction of the system. If we now substitute these 
variables in expressions (|C8|) and (|C10[> . and then in 
Eq. (|6.3p . we arrive at the following relations: 



\Jn 2 - 4 
XA 



(*7„ + XaU), 

= (*/„ + XaIu) - Vn 2 - 4/„, 



(6.5) 



where 



4a - 2 V o 71-0 

XA — e m 2— 



(6.6) 



only depends on the homogeneous variables. Notice also 
that, at the perturbative order of our approximation, the 
substitution in the above relations of the homogeneous 
variables a, 7r Q , tp, and tt v by their barred counterparts 
appearing in Eq. (|3.23p is totally irrelevant. 

The change (|6.5p can be regarded as a Bogoli- 
ubov transformation once we consider the creation and 
annihilation-like variables associated with each pair of 
variables, given by Eq. (|4.1I) and 



6* = 



1 



(6.7) 



together with the complex conjugate. A simple calcula- 
tion allows us to compute the coefficients associated with 
the antilinear part of this Bogoliubov transformation: 



Xa + 3 



2 y/n 2 



lyfn 2 



In order to determine whether the transformation admits 
or not a unitary implementation, the necessary and suf- 
ficient condition is just that 



^T|/3„| 2 =]Td n |/3 n | 2 <cx), 



(6.9) 



nlm 



where d n is a degeneracy factor accounting for the dimen- 
sion of each eigenspace of the Laplace-Beltrami operator. 
Explicitly, we have on the three-sphere that d n = n 2 . As 
a result, the canonical transformation between our vari- 
ables /„ and 7f/ re , and the gauge-invariant ones, \l/ n and 
n , turns out to be implementable as a unitary one at 
the quantum level, since the coefficients j3 n arc clearly 
square summable, including their degeneracy. Therefore, 
the two considered physical descriptions are completely 
equivalent. 



B. Gauge a n = b n = 

Let us now consider the gauge fixing carried out in 
Sec. El Again, one can find the relation in the reduced 
system between the canonical pair /„ and Kf n and the 
gauge invariants ^> n and ILj n . For this, one can first ob- 
tain the expressions of £ ™ l and in terms of /„ and 
nf n (in the gauge under consideration) and then use 
Eq. (|6.3p . A simple calculation yields an expression sim- 
ilar to Eq. (|6.5p . but substituting \A by 



XB 



7i> e 2a 



(6.10) 



We can then compute the corresponding /^-coefficients of 
the Bogoliubov transformation relating the creation and 
annihilation-like variables associated with the variables 
/„ and nf n , on the one hand, and the gauge-invariant 
variables, on the other hand. The result is the same as 
in Eq. (|6.8p . but now with \b replacing \A- Following 
the arguments of Sec. IVI A[ we can easily check the uni- 
tary implementability of this Bogoliubov transformation, 
since condition (|6.9p holds. 



VII. SUMMARY AND CONCLUSIONS 

In this article, we have presented the quantization of 
an inflationary universe with small inhomogeneities. The 
matter content is described by a minimally coupled, mas- 
sive scalar field. First, we have studied the corresponding 
symmetry-reduced homogeneous and isotropic system in 
Sec. |TTJ In order to employ a polymeric quantization Q 
for the geometry degrees of freedom, its classical phase 
space has been parametrized by a densitized triad and an 
Ashtekar-Barbero connection. As for the (homogeneous) 
matter content, we have applied a standard Schrodinger 
quantization. In this kinematical arena, we have intro- 
duced an operator representation for the Hamiltonian 
constraint. The geometric part of this quantum con- 
straint is a difference operator that only relates states 
with support in semilattices with points separated by 
constant steps of four units. Additionally, these semilat- 
tices can be characterized by a continuous, non-vanishing 
parameter e £ (0,4], which is proportional to the mini- 
mum physical volume allowed in the quantum theory on 
that semilattice. Moreover, we have argued that every 
sector of the kinematical Hilbert space with support on 
any of those semilattices is superselected. As any phys- 
ical state belongs to the kernel of the constraint, one 
can see that the restriction to a particular superselection 
sector allows us to entirely determine any solution if its 
initial data on the minimum volume section are provided. 
The final step to quantize the system to completion, i.e., 
to construct the physical Hilbert space, is to equip the 
space of solutions with an adequate inner product. Based 
on the identification of the space of solutions with the 
space of initial data for minimum volume, we have picked 
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up this inner product by selecting a complete set of ob- 
servables and imposing reality conditions on them [30| . 
namely, by requiring that the complex conjugation re- 
lations in this set become adjointness relations between 
the operators which represent the observablcs. Our pro- 
posal allows us not only to formalize the quantization 
of the system, but also to proceed in the analysis of its 
quantum dynamics. The extension of the systematic pro- 
cedures commonly used in LQC to systems without an 
explicit separation of the matter content from the geome- 
try is not obvious, yet our proposal sheds new insights to 
confront also those more general situations. Our prelimi- 
nary numerical analysis shows that, with an appropriate 
choice of initial data, it is possible to recover physical 
states that represent an expanding universe preceded by 
a contracting one with a well-defined quantum bounce 
connecting both of them, even for a massive scalar field. 
Those results are very promising and a more careful study 
will be the subject of future research. 

On the other hand, to progress in the applicability 
of LQC to the analysis of the Early Universe and con- 
struct a quantum framework to study the development 
of small inhomogeneities in inflationary scenarios, we 
have introduced local degrees of freedom in our model by 
means of perturbation theory around the homogeneous 
and isotropic solutions. We have concentrated our efforts 
in the understanding of scalar perturbations, because of 
the fundamental role that they play in present observa- 
tional cosmology. Additionally, from the technical point 
of view, scalar modes are more involved than vector or 
tensor ones, since they incorporate both physical and 
gauge degrees of freedom and their discrimination is not 
trivial. In order to distinguish the true physical degrees of 
freedom, we have carried out two different partial gauge 
fixings. The first one is the so-called longitudinal gauge, 
commonly used in standard cosmology, and the second 
one is the natural gauge fixing in which the spatial met- 
ric is purely homogeneous (i.e., all the inhomogeneities of 
the metric are encoded in the perturbations of the lapse 
and shift functions). In both cases, we reach a descrip- 
tion of the inhomogeneities of the corresponding reduced 
system in terms of the matter perturbations. We have 
introduced a canonical transformation in each case, with 
an eye on the standard Fock quantization of the inho- 
mogeneities, quantization that we have carried out later 
on. Such a new choice of canonical fields is motivated by 
the recent uniqueness results regarding the quantization 
of linear fields (generically) in non-stationary spacetimes 
[25j], where the choice of a Fock quantum theory with a 
vacuum state invariant under the spatial isometries and 
unitary dynamics selects privileged canonical field vari- 
ables for the description of the system, together with a 
specific quantum representation for such field variables. 
We have then adopted this Fock representation for the 
inhomogeneities and combined it with the polymeric de- 
scription initially introduced in Sec. Ill Bl 

Furthermore, we have presented a quantization pre- 
scription for the Hamiltonian constraint in Sec. lIVB] The 



corresponding quantum operator respects the superselec- 
tion sectors of the homogeneous setting. Therefore, any 
state which is a solution to the constraint has support 
in semilattices of constant step in the physical volume, 
starting from a non-vanishing minimum value of it. We 
have been able to prove that, assuming that the solutions 
to the constraint admit an expansion in terms of the per- 
turbative parameter e, the lowest order contribution in 
the perturbative expansion of a solution can be regarded 
as a background state, in which the inhomogeneities play 
no dynamical role, and which can be totally determined 
by its initial data on the volume section v = e. On the 
other hand, it can be seen that the next contribution 
in the expansion of the solution satisfies the FKW con- 
straint equation but with a source term coming from the 
interaction of the inhomogeneities and the background 
state. Again, we have proven that this contribution to 
the solution can be determined once its initial data on 
the minimum volume section are provided. Therefore, 
in order to construct the physical Hilbert space, we just 
have to endow the space of initial data for the solutions 
with a Hilbert space structure, supplying it with an inner 
product which, for instance, can be selected by requiring 
reality conditions on a complete set of observables de- 
fined on the section v = z. The very same procedure has 
been carried out in Sec. |V]for the alternate gauge fixing 
considered in this article. 

In addition, to avoid any gauge-fixing dependence of 
the possible physical outcomes of our proposals, we have 
introduced a family of gauge-invariant variables. We have 
constructed a canonical pair of conjugate variables out 
of the gauge-invariant energy and velocity perturbations 
(see Ref. [21| and Appendix [C]). Adopting a standard 
Fock description for them (namely, the so-called mass- 
less representation), we have been able to prove that the 
fundamental matter perturbations in each (gauge-fixed) 
reduced system and the new gauge- invariant variables are 
related by means of a canonical transformation that can 
be implemented unitarily quantum mechanically. More- 
over, the uniqueness results of Ref. pBj can be immedi- 
ately applied to the gauge-invariant variables proposed 
in this article. Therefore, the physics predicted by the 
different descriptions proposed in this work (either based 
on gauge invariants or not) is equivalent, at least as far 
as standard quantum field theory in curved spacetimes is 
concerned. 

In conclusion, we have been able to provide a full quan- 
tum description of an inflationary universe with small in- 
homogeneities propagating on it, in the context of LQC. 
The model is now ready to produce physical predictions, 
which will be the aim of future work. 
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Appendix A: Unitary transformation 

In this Appendix, we will provide a unitary transfor- 
mation that cancels the complex phases in the definition 
of the operator f2 2 constructed from Eq. (|2.4[) . This op- 
erator is given by 

+ (\R + (v)\ 2 + \R-(v)\ 2 ), 



(Al) 



where 



R±(v) = ^^|«| 1/2 |« ± 2I 1 / 2 [sgn(«) + sgn(v ± 2)] 



xe 



(A2) 

and p,(v) = [A/^Ghv)} 1 / 3 . 

Given a generic function h(v), we define its complex 
exponentiation 



(A3) 



We consider now the operator Uil 2 U 1 and, in particu- 
lar, the N^n term (the contribution of N-4,n is given by 
the adjoint, and the remaining part of ft 2 is a multiplica- 
tive operator in the v- representation) . We can see that 
the complex phases cancel if, Vw, 

h{v + 2) - h(v - 2) = £<E+2 + + m . (A4) 

Let us change the label v — > v + 2 and restrict the trans- 
formation to a superselection sector with, e.g., positive 
orientation of the triad, i.e., v — e + 4n (n € N), without 
loss of generality. Then, we are able to determine the 
explicit form of the function h(y) by a recursive process: 



h{e) 



H{s) 



ft( e + 4) = ME+i) +A(e + 2) + m + h(e ) 



fi(e + 4) 



+ /2(e + 2)+/i(£), 



, / on m( £ + 8) _, „. fi(e + A) A . 
h (.£ + 8) = — - + (i(e + 6) + ^— — - + h(e + 4) 



p.(e + : 



+ £(e + 6) + /2(s + 4) + ^(g + 2) + /2(e), 



ft(g + 4n) = 



H{e + in) 



2n— 1 



^ /2(e + 2j), n e N+. (A5) 

3=0 



The constraint C, involving the operator ft 2 , superse- 
lects sectors. Once the analysis is restricted to a specific 
sector, it is natural to introduce in it a unitary transfor- 
mation of the above form. In this sector, the operator Q 2 
is then mapped into f2§ [see Eq. (|2 .6[) ] . It is worth com- 
menting that, for large n, the function h(v) converges to 
the function v 2 ^ 3 (up to a factor), essentially because the 
sum in (|A5I) converges to the integral J dv v^ 1 / 3 ~ u 2 / 3 
and the first term in that expression is fl(v) ~ o(v — 1 ' 3 ). 
Therefore, we recover the function of the corresponding 
unitary map introduced in Ref. Q. 



Appendix B: (Hyper-)Spherical harmonics 

We now briefly summarize the main properties of the 
(hyper-)spherical harmonics Qf m on S 3 [3ll l36l|. They 
form a basis (of square integrable functions with respect 
to the volume element defined by the standard metric on 
the three-sphere) in which the Laplace-Beltrami opera- 
tor is diagonal, with a discrete and unbounded negative 
spectrum. They are labeled with three integers: n, I, and 
to. The last two account for the degeneracy of each eigen- 
value of the Laplace-Beltrami operator. Their ranges are 
< I < n — 1 and —l<m<l. We will display these la- 
bels explicitly only in those steps of our analysis in which 
they play a relevant role. 

Thus, the scalar harmonics Q n are eigenfunctions of 
the Laplace-Beltrami operator on S 3 , and specifically 
they satisfy 



\a 



\a _ 



(n 2 -l)Q n , n= 1,2,3. 



(Bl) 



Here, the symbol | denotes the covariant derivative with 
respect to the round metric £l a b on the three-sphere of 
unit radius. Since we deal with real scalar fields in our 
discussion, we choose the harmonics Q n to be real. 

Starting from these scalar harmonics, it is straightfor- 
ward to construct a family of vector harmonics P™ by 
applying covariant derivatives: 



pn _ 



1 



1 



2,3,4. 



These vector harmonics satisfy 

> 2 -3)p;, p«i° 



pn \b 
r a \b 



(B2) 



-Q n . (B3) 



In addition, we can also construct a family of tensor 
harmonics, namely, the scalar tensors 



Qab = l^abQ n , n= 1,2,3. 



(B4) 



and the traceless tensors 



P ab = -^Q n \ab + \n ab Q n , n = 2,3,4... (B5) 
n A — 1 1 3 
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These traceless tensors have the following properties 



(n 2 -7)P™ 6 



c 6 |b --3(« 2 -4)pr, 

P? b lab = l(n 2 -4)Q n . 



(B6) 



Finally, if we call dv the integration measure on S 3 
corresponding to the volume element determined by the 
metric f2 a &, and we normalize the scalar harmonics so 
that 



/ 



d v QfmQvm' — Snn'dll'Sr, 



(B7) 



it is straightforward to check that 



dv(P a )UP a )?'m> 



n 2 - 1 



3nn'&ll' $mm' ? 



dv(Q ab )UQ ab )u m > = ^ nn .5 u ,5 mmli 

, / 2 n 2 — 4 
dv(P ab )UP ab )?<m' = lnr -rS nn ,5 u >6 mm ,. (B8) 
3 n z — 1 



Appendix C: Bardeen potentials 

In this appendix, we provide the definitions of some 
relevant gauge-invariant quantities [2lj and discuss the 
dynamical relations between them. In doing this, we use 
the expansion of the metric components and the matter 
field given in Eqs. (I3.1[) - (|3.4[) . We consider gauge trans- 
formations of the form x'^ L = + e£ M , introduced in 
Sec. IVI1 with the parametrization (|6.1[) and (|6.2p . 

After a transformation of this kind, the modes of the 
metric can be written in terms of the original ones and 
the modes of 



9n >-» 9n + -^£o, 

N 



(CI) 



i /l 



in «-> a„ + — -C + aCo , 



bn •->■ b n - — f, 
3e Q 

fn ^ /n + ~£o ■ 



(C2) 



Recall that the overdot stands for the derivative with 
respect to the conformal time 77, and u) 2 = n 2 — 1. With 
these transformation rules, the following scalar modes 
define gauge-invariant quantities: 



1 d f e a 



e a ujl dr] V No 



-k n 3e 6 n I , 



* r. 
f" 



a n + b n + ■% ( -77- - 36„ ) , 



(C3) 
(C4) 



1 



E e 2 



<pfn - <p 2 g n + {iaip + e 2a m 2 tp)f n (C5) 



u> n \ if N 



(C6) 



where E — (e 2a tp 2 



rh 2 (p 2 )/2 is proportional to the 



energy density of the background scalar field. If we com- 
pare these quantities with the gauge-invariant variables 
originally defined in Ref. 21], we see that the gauge- 



invariant $^ can be identified with the potential defined 
in Eq. (3.9) of that reference. On the other hand, the 
Bardeen potential defined in Eq. (3.10) of that work cor- 
responds to the quantity $ s . Our invariant energy den- 
sity perturbation, £ m , is the gauge invariant of Bardeen's 
Eq. (3.13). Finally, the gauge- invariant matter velocity 
of Eq. (3.11) in Ref. [2l[ corresponds to v s . 

Let us derive now the expressions of £™ and v n as func- 
tions of the canonical variables introduced in Sec. IIIIl For 
the gauge-invariant energy density and velocity pertur- 
bations, we need to employ the dynamical equations in 
order to rewrite the momenta 7r/„ and ~K bn in terms of the 
time derivative of the corresponding conjugate variables, 
i.e., f n and b n , respectively. For the matter perturbation 
we obtain 

?r/„ = e 2a f n + 7r v (3a„ - g n ). (C7) 
Taking this into account, a simple computation yields 
1 



On the other hand, for the momentum conjugate to b n 
we have 



1 



> 2a L 



inabn ~r 2a k)- (C9) 



Therefore, the velocity gauge-invariant perturbation can 
be written as 



1 



" 2 f " 

UnJn 2a~ 

V 



n 2 - 1 



7T bn + 4-Kabn 



;t < -" \ n — 4 

'(CIO) 

We now turn to the relations existing between these 
gauge-invariant quantities. On the one hand, one can 
check that the equation of motion for the perturbation 
b n (see Eq. (Bll) of Ref. [3l|) is equivalent to 



^ + ^ = 0, n = 3,4. 



(Cll) 



This equation corresponds to Eq. (4.4) in Ref. [21], valid 
for any fluid whose stress-energy tensor has a vanishing 
traceless part. Another interesting relation arises from a 
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linear combination of the constraints if," and fFj , given 
in Eqs. f|3 . 10[) and (|3.11[) . respectively. One can see that 



H![ - iaW\ - 3H lo a n = e ia E £, 



1 



e a (n 2 -4)&i 



If the conditions HT\ = and = are satisfied, 



the last equation reduces to Eq. (4.3) of Ref. [21 
equivalently, 



or, 



' n 



Let us now define 



Pn 



3e 



2<i 



771 cm 

;^o£ n 



Vn 



1/V 

2 \e 2a 


~ 2 2\ 

- m tp j, 


Po 
Eo 


dP 


2e 2a rh 2 

1 + ^7-^- 

daip 




dE 




SPn 

Pa 


dP SE n 
dE Po 


1 -c 2 

x l -s cm 
c n 

W 



(C12) 

(C13) 
(C14) 
(C15) 



Here, Po is proportional to the pressure of the homoge- 
neous matter field, while 5E n and 8P n are the modes 
of the energy-density and pressure perturbations. Us- 
ing these formulas, one straightforwardly proves that 
Eq. (4.5) of Ref. [21] can be written as 



rp cm jr 

h ° b n - ^ 



--f( e X) 



e a <S>C 



(C16) 



Finally, we can obtain an energy equation like Eq. (4.8) 
of Ref. 21]. We start from the equation of motion of /„, 
which is a second-order differential equation that can be 
found in Eq. (B14) of Ref. [3l|. If we combine it with 
Eqs. (|C11[) and (IC12I) . and with the equation of motion 
of a, we arrive at the expression 



d_ 

d:q 



1 e 3a En£L 



-e a (p 2 uj n v s n 



(C17) 



Appendix D: Gauge-invariant curvature 
perturbation 

Another interesting gauge-invariant quantity is the 
curvature perturbation, which is defined as 



Tin = - — < 



(Dl) 



and was originally studied by Bardeen (see the corre- 
sponding definition in Eq. (5.19) of Ref. |2l|). In the case 
of spatially flat cosmologies, it is common to work with a 
closely related quantity, the so-called Mukhanov- Sasaki 
variable v n = zTZ n , with z — e a tp/a. In such situations, 
the power spectrum of primordial perturbations can be 
easily derived using the Mukhanov-Sasaki variable. 

In our case, this gauge invariant has different expres- 
sions in terms of the canonical variables employed in the 
two distinct gauge fixings considered here. 

a. Longitudinal gauge 

In order to compute 1Z n in terms of our fundamental 
fields f n and w/ n , we recall first that the potential $^ is 
related with £™ by means of Eq. (|C12j) . Applying the 
results of Sec. IVI| it is then easy to obtain that 



37T„ 



e 3a (n 2 - 4) 



(*/„ + XA.fn) + 



(D2) 



The variable \A was defined in Eq. (|6 .6[) . Clearly, in 
the large n limit, the main contribution to the previous 
expression comes from the last factor, proportional to f n - 



b. Gauge fixing a„ = b n = 

This gauge provides a simple relation between the cur- 
vature perturbation and the modes f n . Actually, using 
the definitions (|D1|) . (|C4[) , and (|C6I) . a straightforward 
computation yields 



T^-n — fn- 
e a 7r M 



(D3) 



Therefore, in this gauge, where the spatial geometry is 
homogeneous, the unique contribution to the curvature 
perturbation comes essentially from the perturbation of 
the scalar field. 

Moreover, in the ultraviolet limit, the form of the cur- 
vature perturbation coincides formally with that of the 
longitudinal gauge, given by Eq. (|D2I) . As a consequence, 
for the two studied gauges, the Mukhanov-Sasaki vari- 
able v n either coincides with /„ (gauge a n = b n = 0) or 
converges to it (longitudinal gauge). 
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